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The large N sigma model, in D < 4 space-time dimensions, with disorder 
a function of d space dimensions, is analyzed via a renormalization group treatment. 
Critical exponents for average quantities are calculated, first to lowest order and then to 
all orders, m e = D — 2 — |. In particular, it is found that vd — 2. When D = d+ 1, this 
model is equivalent to a large N limit of the strongly commensurate dirty boson problem. 

I. INTRODUCTION 

The problem of the statistical mechanics of a system with quenched randomness is very interesting 
and difficult, combining the problems of disorder and interaction. We can consider two cases, the ordinary 
statistical mechanics of such a system, or the quantum statistical mechanics, in which a d dimensional 
quantum system can be considered as a d + 1 dimensional classical system, with disorder correlated in 
the time direction. These correlations make the effects of disorder much stronger, as seen, for example, 
in the one dimensional disordered quantum Ising chain [jj], a problem which is probably the most well 
understood disordered quantum phase transition. There has also been much work on the quantum 
statistical mechanics of a system of repulsively interacting bosons in a disordered potential, known as the 
dirty boson problem , but this problem is not as well understood. 

An action that may be used to describe the dirty boson problem is 

J d d x dt (-da${x, t)d x <j>(x, t) - d t ~$(x, t)^(x, t) + w(x)^(x, t)d t <j)(x, t) - U{x)M> + jj(W) 2 ) W 

Here, U(x) and w(x) are considered to be quenched random variables, with small fluctuations. Although 
the action above describes a D — d+l dimensional system, the randomness in U(x) and w(x) is a function 
of the d dimensions x only. This action is expected to describe a system with a number of different phases. 
If the average chemical potential U is negative, the system is a Mott insulator, which has a gap. Increasing 
the chemical potential, the system then becomes a Bose glass, and finally a superfluid. The Bose glass 
is characterized by a gapless phase with infinite superfluid susceptibility fic|| . In this phase, the bosons 
occupy localized states, with a localization length that diverges as the superfluid transition is approached. 

If the boson density commensurates with a lattice, then w(x) vanishes on average; this can 
be accomplished experimentally by tuning a chemical potential. In a typical experimental situation, 
however, there will still be fluctuations in w(x). This case, with nonvanishing fluctuations in w(x), is 



the weakly commensurate case. If w(x) = identically, then we have the strongly commensurate case. 
In an experimental setup with either a dirty boson system or a Josephson junction array, this would 
require tuning the local chemical potential to be constant, so there is no local breaking of particle- hole 
symmetry due to disorder, while still having some fluctuations in the hopping term. This is obviously a 
more difficult task. 

In this paper, we consider only the strongly commensurate case, where w(x) = everywhere. 
Although less physically applicable to boson systems, we feel that the results on the strongly commen- 
surate case are interesting in themselves, as well as being applicable to some other problems of quantum 
critical points. In particular, there exist quantum critical systems, with an 0(3) order parameter, which 
in the disordered case may be well described by an N = 3 disordered rotor model like those consider in 
this paper. See work by Sachdev [Q, and references therein to collinear quantum antiferromagnets. 

In the strongly commensurate case, we will not see the Bose glass phase. As the chemical 
potential is increased, the system will go from a gapped Mott insulator phase to a gapless Griffiths phase, 
but the superffuid susceptibility will remain finite. One reason for considering the strongly commensurate 
case is related to the use of the large N limit and will be discussed later. Results will be described in a 
future work for the weakly commensurate case j|. However, we feel that it is useful, in this paper, to 
mention the existence of the weakly commensurate problem, and the Bose glass found in that problem, 
to understand the different phase diagram with and without the w(x) term. The Bose glass phase, as we 
understand it, relies on the existence of the w(x) term, and so, contrary to what other authors suggest 
we feel it is impossible to have such a phase in the strongly commensurate problem. Since the problem 
without the w(x) term can be understood simply as a rotor problem, it seems that the only possible 
phases are paramagnetic and ferromagnetic, with a Griffiths phase near the critical point. 

In this paper, we consider a large N generalization of the strongly commensurate dirty boson 
problem [0. We consider a system described by the following partition function 



where the field 4>i has N components, and is a function of d space coordinates x and D — d time coordinates 
t. The (5-function constraint on the length of the spins is technically simpler and is not expected to alter 
the universality class from the quartic interaction considered above. The function o 2 {x) is a function 
of the spatial coordinate x only, and is related to the spatially varying U (x) considered in the previous 
equation. 

The advantage of the large N formulation of the problem is that the system may be exactly 
solved for any fixed realization of disorder by solving a self-consistency equation. We may replace the 
(5-function constraint by an integral over a Lagrange multiplier field X(x) as 
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where the integral for A(x) extends from — zoo to +ioo. In the large N limit, we can use a saddle point 
approximation for A, the saddle point being found by the self-consistency equation 

a 2 (x) = (x,t = 0\(-dl + X(x))- 1 \x,t = 0) (4) 

This equals 

a 2 (x) = J d {D - d) oj(x,t = Q\(-dl + oj 2 + \(x)) \x,t = 0) (5) 

We will write a 2 = a 2 + 5a 2 , where Cq is a constant, which is tuned to drive the system through the 
phase transition, and Sa 2 is a random term. For small a 2 the system is in the Mott insulator phase. For 
large CTq the system is in the superfluid phase. 

In the above equations, we assume that the Green's function on the right-hand side has been 
renormalized by subtracting a divergent quantity. That is, we will take a Pauli-Villars regularization for 
the Green's function, and take the regulator mass to be very large, while adding an appropriate divergent 
constant to a 2 on the left-hand side. The cutoff for the regulator is completely different from the cutoff 
for fluctuations in 5a 2 that will be introduced for the RG of the next section; the cutoff for the regulator 
will be much larger than the cutoff for fluctuations in da 2 and will be unimportant in the RG. 

Note that in the large N limit, it makes sense to define the theory for non-integer D — d, by 
doing the integral over u>. For finite N, such a definition may have trouble, and indeed in the double 
dimensional expansion in 4 — d and D — d, there is some question about expanding around D — d = § . 

One disadvantage of the large N expansion is that it is slightly difficult to include the terms 
linear in the time derivative, necessary to understand the weakly commensurate and incommensurate 
dirty boson problem. Consider the simple problem 

where S is a small number introduced for point-splitting. While for U < everything is correct, for U > 
the solution of the problem via the self-consistency equation becomes ill-defined for zero temperature, 
although not for arbitrarily small, non-zero, temperatures. The trouble is in the zero temperature limit of 
the problem in which we must consider frequencies u> arbitrarily close to 0. This will be further discussed 
in future work B. 

Having dropped terms linear in the time derivative, we do not expect to see a Bose glass phase. 
The original argument for the Bose glass phase was based on considering a system consisting of localized 
states for the bosons, with hopping between the localized states being neglected. It was then shown 
that for a random distribution of chemical potentials for each localized state, the spectrum of excitations 
would be gapless with constant density of excitations of low energy. These excitations would correspond 
to changing by one the number of particles in a given localized state. However, if one considers a system 
of localized states with action containing terms quadratic in the time derivative instead of linear, in the 
large N limit, the density of low energy excitations goes to zero. In the Griffiths phase of the model we 
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are considering there are states at arbitrarily low energies, but the density of states vanishes as e~ , 
where c is some constant. 

In the RG treatment, we will be considering the phase transition between the Griffiths phase 
and the superfluid phase, with weak randomness in a 2 (x). By the Harris criterion [T^ ], weak randomness 
is irrelevant at the pure fixed point for vd > 2 and marginal for vd = 2. There is also a bound that 
for a stable disordered fixed point vd > 2 |l3). For D < 4, we find that v = 1/(D — 2). Thus, there is 
a range of values of D and d satisfying d/(D — 2) = 2, at which the randomness is marginal. We are 
able to construct a renormalization group near any of these values. Note that for the quantum statistical 
mechanics of a d dimensional system, where D = d + 1, we find that the disorder is marginal for d = 2 
and relevant for d > 2. The renormalization group is not constructed by expanding down from an upper 
critical dimension. Instead, it is constructed by expanding upwards near a range of dimensions such that 
dj (D — 2) = 2. The RG is only valid for D < 4, as discussed in the next section. 

Within the RG, we will proceed perturbatively, but will obtain exact results for all the various 
exponents for average quantities. One of the most striking results is that vd = 2. This implies that the 
system saturates the bound discussed above jl3|] . Similar results have been found for various other phase 
transitions. For example, exact results for the transverse field Ising model in d=l lead to saturation of 
this bound |IJ], with v = 2. Also, numerical simulations on the quantum Ising spin glass in 3+1 dimensions 
show this bound This may be a common feature of quantum phase transitions. 

The next two sections will develop the renormalization group. The procedure will be to take a 
problem with fluctuations in a 2 up to some cutoff A, and define another problem which includes fluc- 
tuations in a 2 only up to a wavevector A — <5A, such that we preserve the low-momentum correlation 
functions. In section II we will perform a one-loop Wilson-Fisher RG. To extend this to higher orders 
would require keeping track of many operators, while an alternative procedure discussed section III re- 
quires only keeping track of the renormalization of two quantities: a 2 and S, where S is a measure of the 
strength of fluctuations in the disorder Sa 2 . In order to preserve the low- momentum correlation functions, 
we will require that the propagator used includes self-energy corrections due to the high- wavevector fluc- 
tuations, and also require a renormalization of the disorder strength to produce the same low-wavevector 
fluctuations in A, up to a renormalization of the vertex. Requiring the same low-wavevector fluctuations 
in A makes the procedure consistent, and permits one to continue iterating the RG. 

II. RENORMALIZATION GROUP 

We proceed with a renormalization group acting directly on the self-consistency equation. The 
goal will be to start with a given problem which includes high wavevector fluctuations in a 2 , and find a 
related problem which includes fluctuations in a 2 only up to a lower wavevector. We will require that 
the correlation functions in the related problem are equal, up to rescaling, to the correlation functions 
in the original problem averaged over the high wavevector fluctuations. This requirement is what will 
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define the RG. 

At lowest order, considered in this section, this will require that A is unchanged, up to rescaling, 
at low wavevectors. To this order, we can accomplish this goal by defining a new problem in the same 
form as the old problem, via a self-consistency equation with one parameter A; to higher orders this 
Wilson-Fisher procedure will be more complicated and we will instead use an alternative technique in 
the next section. To higher orders within the Wilson-Fisher approach one encounters operators such a 
spatially fluctuating gradient terms and to 2 terms, as well as non-Gaussian and momentum dependent 
distributions of the disorder. 

In equation (^J), we have broken a 2 (x) into two pieces: a 2 + 5a 2 (x), where Sa 2 (x) has vanishing 
mean. In the perturbative e expansion being developed here, we can perturb in 6a 2 (x) as it will have 
fluctuations at the critical point which are of order e. 

Consider the self-consistency equation (^) , and define a cutoff A such a 2 only has fluctuations for 
wavevectors less than A. Formally, we can invert equation to obtain A as a function of a 2 as follows. 
For small 6a 2 , we can expand the right-hand side of (|j) as a power series in A. At criticality, where A 
vanishes for vanishing Sa 2 , we find to lowest order 

Sa 2 ( P )= [ d d kd D - d cj- -i -_L_A(p) (7) 

K ' J (p - k) 2 +lo 2 k 2 + to 2 y ' v ; 

This can be derived from a diagram similar to a polarization bubble as shown in figure 1. 
This implies that as a formal inverse, to lowest order, 

X(p)=c 1 5a 2 (p)p i - D (8) 

where c\ = r(2-_D/2)r li (D/2-i) • ^ or smau wavevector components of Sa 2 , this formal inverse is ill- 
behaved, but for high wavevector components when D < 4, the procedure given will be correct to lowest 
order in e. This is why we need D < 4 for the RG, as otherwise the polarization bubble is divergent for 
large wavevector instead of small wavevector. 

Define a measure S of the strength of disorder, by assuming that 

(Sa 2 (p)6a 2 (q) = (2Tr) d 6(p-q)S (9) 

where to lowest order we may take 5a 2 to be Gaussian distributed. Then, we find that at the cutoff A 
has mean-square cfA 8 ~ 2D (2n) d S . We will define L = c\A 8 ~ 2D S to measure fluctuations in A. The effect 
of these components in A will be to renormalize the propagator for the field cf>, as well as to renormalize 
the vertex used to calculate the scattering of off A and to calculate a 2 . The self-energy is given by 

Z( P ,u) = 6a[ d^k 1 2 L (10) 

See figure 2 for the appropriate diagram. This is equal to a constant, which may be absorbed into a 
renormalization of A, plus 
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S -±(c 2P 2 + cW)L + ... (11) 



where C2 = (1 — 4/d)c3 and c 3 = 2 r(d/2) ^ 4 ' This implies that the propagator is renormalized to 

1 



(12) 



(1 + fc 2 L)p 2 + (1 + f C3 L) W 2 
while the vertex renormalization, shown in figure 3, is given by 

— c 3 L (13) 

This implies that we can write a new self-consistency equation 

(1 - S -±c 3 L)a 2 (x) = J S d -^uj (x,t= 0| ( (1 + 5 -±c 2 L)dl + (1 + ^c 3 L)(c 2 + \{x))) ^)\x, t = 0) 

(14) 

where a 2 now has only components at wavevectors less than A — SA. Note that by removing the high 
wavevector components of A in the vertex renormalization, we reduce the value of a 2 on the left-hand 
side of equation (|l4|). We rescale equation ( |l4| ) by defining a new scaled A, scaling the integral over u>, 
and scaling the d space dimensions, to obtain 

(1 + 5 -£(D - 2))(1 - ^(c 3 - c^Lf-^^ix) = J A D -*>w (x,t= 0\(-0 2 x + oj 2 + X(x)y 1 )\x,t = 0) 

(15) 

Assuming that a = ctq + 5a 2 , with Sa 2 having a distribution given by equation (|^), we can obtain 
RG equations for a 2 and S by considering the scaling of a under equation (Eq). We obtain 



^_ B - 2 - to -«tfA«»5(I-^) (16) 

\^- D -'-^-'^- 2D ^- £ ¥>-' 1 ' 2 < 17 » 

Therefore, there is a fixed point of the RG flow, for D - 2 - d/2 = (c 3 - c 2 )c? A 8 - 2D 5(1 - ^). 
The fixed point is stable in the S direction. It is unstable in the ctq direction. We find 9 — d/2 at the 
fixed point, leading to the exponent v — 2/d, which implies that the bound vd > 2 is saturated. 

In fact, equation ( |l6|) is not quite correct. Due to scaling of the regulator mass, there is an 
additional constant term in the change of a 2 with respect to InA, or equivalently an extra term in the 

CT 

change of lnc 2 , with respect to InA. Here, C is some non-universal constant. This additional term only 
leads to a non- universal change in the critical value of <Tq, and no change in the critical exponents. 

We can obtain the dynamic critical exponent by considering the different scaling of p 2 and uj 2 in 
the propagator. This difference is (c 3 — C2)c 2 A s ~ 2D S = ° j ■ Therefore, 



G 



z = l + 



2-D + d 



2 + d-D 



d/2 



(18) 



The average Green's function, (G(p,uj)), is 



1 



(19) 



P 



2 + d~D _|_ w 2 



This result for the dynamic critical exponent differs greatly from that expected in the weakly 
commensurate case, where scaling theory predicts z — d This is not due to 1/N corrections, but 
rather a big difference between weakly and strongly commensurate phase transitions j3[. 



Given the simplicity of the result vd = 2, one might suspect that the exponents for average 
quantities remain unchanged to all orders. A very simple argument shows that this is true. To obtain 
results to higher order, it is useful to use a formulation of the RG other than the Wilson-Fisher RG. 
We will proceed as follows. Consider a theory with quadratic fluctuations in a 1 of some magnitude S, 
with fluctuations ranging up to some wavevector A, and with given ctq. In this theory calculate the 
average Green's function at low momentum, as well as the low momentum fluctuations in A, to some 
order in perturbation theory. The reason for considering fluctuations in A is that this is what is needed 
to continue the computation of the Green's function to higher orders. Then, for a theory with cutoff 
A — SA, determine the appropriate value of S and appropriate coefficients of d% and 9 t 2 in the action, so 
as to reproduce the given low energy behavior to the same order in perturbation theory, with the same 
A up to a renormalization of the scattering vertex. 

This procedure is a well-established alternative to the Wilson-Fisher and Callan-Szymanzik tech- 
niques. It is discussed for example in the classic reference Domb and Green volume 6, as the second of 
two field theory techniques for performing an RG Q|. Within this procedure, one needs only to keep 
two terms in the RG, the value of (Tq and the value of S, as these are the only two relevant terms. It is 
possible, of course, that higher loop corrections will make other operators relevant, and destabilize the 
fixed point; this is something that cannot be analyzed within this approach. However, this approach 
will yield results that agree with the Wilson-Fisher technique so long as the fixed point is stable. Let 
me emphasize this point: if no other operators become relevant, then this technique agrees with the 
Wilson-Fisher procedure to all orders. If other operators become relevant, then the procedure does not 
work and the fixed point has more than one unstable direction. What we show in this section is that if 
no other operators become relevant and the fixed point remains stable then the exponents are unchanged 
to all orders. What was shown in the previous section is that there is a stable fixed point to lowest order 
for D — d < 2, with exponents as given above. It is reasonable to assume that the fixed point remains 
stable to all orders, at least for some range of e, so that the exponents remain unchanged to all orders at 
this fixed point. 



III. RESULTS TO ALL ORDERS 
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In the case of the one loop calculation above, it was found that the vertex renormalization for A, 
the vertex renormalization for a 2 , and the renormalization of uj 2 all had the same coefficient. This will 
remain true to all orders. This is due to a Ward identity, explained at the end of this section. Therefore, 
the general RG equation to all orders will take the form 

(1 - S ~^F 1 (S, D, d)))a 2 (x) = J d^-^u (x, t = 0| (-(1 + 6 -^F 2 (S, D, d))8 2 x + (1 + ^F^S, D, d))(co 2 + A(.t))) ^)\x, t = 0) 

(20) 

where F\ , F 2 are some generic functions of disorder strength and dimensionality. By the same procedure 
of rescaling, and writing a 2 = o\ + 5a 2 , as before, we can obtain 

^. = D-2-(F 1 (S, D, d) F 2 {S, D, d))(l - ^) (21) 

l^x = D 2 {Fl {s > A d) F2{s > A d)){l ~ d/2 (22) 

Then, although the fixed point of the RG flow will be at a different value of S than was obtained to 
lowest order, the exponents will be unchanged from their lowest order values. 

Let us consider why this Ward identity holds. First, consider the equivalence between the renor- 
malization of g 2 on the left-hand side of the self-consistency equation and the renormalization of A on 
the right-hand side. It is apparent that both rcnormalizations are vertex renormalizations and are the 
same by definition. 

The more interesting aspect of the Ward identity is the equivalence between the renormalization 
of a 2 and the renormalization of the lo 2 term. This, however, is again almost true by definition. The 
term lo 2 is constant in space, and a function only of frequency. The Green's function for given frequency 
uj is constant under a shift uj 2 — > uj 2 + A and A — > A — A, for some constant A. For this property to 
hold under renormalization, we need the desired Ward identity. Alternatively we may say that for any 
given frequency ui, the term uj 2 plays exactly the same role in the Green's function as a constant term in 
A would, and so the renormalizations must be equal. 

For a more diagrammatic derivation of the identity, proceed as follows. The coefficient of the lo 2 
term can be obtained by differentiating the inverse of the average Green's function with respect to lo 2 at 
u> = 0. In any given diagram for the average Green's function, one can differentiate the inverse of any 
one of the propagators in the diagram with respect to lo 2 . This gives one a sum over different places to 
insert the derivative. However, each place one inserts gives a result for the diagram which is equivalent 
to placing a scattering vertex for A at that point, and thus the total renormalization of the uj 2 vertex is 
the same as the renormalization of the A vertex. For example, in figure 2 for the one-loop self-energy, the 
renormalization of the uj 2 term can be obtained by differentiating the inverse of the propagator in the 
loop with respect to to 2 . However, this yields exactly the diagram in figure 3. 

The existence of this Ward identity docs not require the formulation of the RG used in this 
section. A similar identity would exist within a Wilson-Fisher renormalization group. This would then 
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connect the renormalization of the to 2 term to that of the a 2 and A terms at lowest order as found in 
the last section. At higher orders, other terms would also be connected. For example, the existence of 
the symmetry would also require that a term such as A(x)ui 2 , that is, a spatially fluctuating lu 2 term, 
would have a renormalization connected to that of a term A(x)X, and a term A(x)<j 2 . This last term 
would imply a spatially varying magnitude of fluctuations in the disorder 8a 2 , which would mean a non- 
Gaussian distribution of the disorder. From the Ward identity, we would still have a connection between 
the renormalization of the lu 2 , a 2 and A terms, which would be enough to obtain the desired results for 
critical exponents, to all orders. 

The only thing the Ward identity does not guarantee is the stability and existence of the fixed 
point. It only guarantees the exponents if the fixed point exists and is stable. As D — d — > 2, the fixed 
point runs to stronger and stronger disorder. This will be discussed more in the conclusion. 



IV. A HEURISTIC TREATMENT OF THE RG 



We present a simple alternative treatment of the problem which supports the results of the 
RG and provides a physical motivation for it. First, let us present an alternative version of the Harris 
criterion. For the pure system, at criticality, the self-consistency equation is given by 

a 2 = effcd^-* (23) 



where we have now inserted a cutoff for high energy states of the field <j>. We find by doing the integral 
in equation ( |23] ) that a 2 oc A^ -2 ). However, if we consider a disordered system, we find that, after 
removing states with k 2 + lu 2 > A 2 , we are considering a case in which we have smeared out space over 
a length of order 1/A, which corresponds to a d-dimensional volume V of order 1/A d . The average value 
of the fluctuation in a 2 over this volume is of order V" 1 / 2 , or A d / 2 . When the Harris criterion indicates 
that disorder is relevant, we find that K d / 2 > K^ D ^ for small A. This is to be interpreted as saying that 
there are not enough low energy states in the pure system to produce the required a 2 . 

However, we expect that disorder will increase the low energy density of states, and we heuris- 
tically modify the Green's function to fca * 2 , with a > 2, and we change the cutoff to k a 
that we have 



a 2 



d d kd^ D ^Lu - 1 „ (24) 
k a + uj 2 y ' 



In this case, we find by doing the integrals that a 2 oc A 2rf /°~ 2 ~ d+ - D . However, we are now smearing the 
system out over a length of order l/A( 2 /°). For small A, this is less the length 1/A which we had in the 
pure system. This is to be expected, indicating the motion has become subdiffusive. This leads to an 
average of a 2 over the d-dimensional volume of order hS d / a ) . Equating these two results for a 2 , one given 
by the integral in equation (124) and the other given by the volume average of a 2 , we find that 
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2d/a-2- d + D = d/a (25) 

or a = d/(2 + d — D). This agrees with the value of a obtained by the RG treatment of the previous 
section. Further, if we imagine slightly leaving the critical point, we can imagine that the self-consistency 
equation gets replaced by a 2 = fca+ ^ +A , where A and a 2 are now taken spatially constant. If we slightly 
adjust a 2 away from the critical value, we can calculate the change in A, and thus the correlation length 
of the system, and we find that we obtain vd = 2. 

V. THE ZERO ENERGY WAVE FUNCTION 

As we increase cTq, with fixed randomness 5a 2 , we find a sequence of different phases. Assume 
we are given X(x) as a function of (Jq + 5a 2 . For small Cq, the operator (—d 2 + X(x)) has a gap, and there 
are no eigenstates of this operator with eigenvalues below the gap. For the case of the zero temperature 
quantum phase transition, where D = d + 1, this gap corresponds to an energy gap for the system. 
This is the Mott insulator phase discussed in the introduction. As (jg is increased the gap decreases and 
at a certain point the system enters a Griffiths phase. In this case, there are eigenstates at arbitrarily 
small energy, all of which are localized, and the average correlation function decays exponentially. At the 
critical a 2 , the average correlation function acquires power law behavior, and there appears a state at 
zero energy which is delocalized. Above the critical (Jq, particles begin to condense into this state. We 
would like to examine the wave function a(x) of this state. 

We will first give a physical argument for the form of the wave function. Then, we will derive 
this to first order in e through the RG. Unfortunately, this result cannot be derived to all orders in e. 

Arguing physically, the following should be a good description of this wave function. Consider a 
problem at short distances. From the RG, we know that, on a short scale, one can obtain A directly from 
a 2 . So, the wave equation, (—d 2 + X)a(x) = 0, can be solved approximately on a short scale knowing 
only the local fluctuations in a 2 . However, the wave equation is homogeneous, and defines a(x) only up 
to a multiplicative constant. So, the zero energy wavefunction in some region in space is defined by the 
local disorder, up to a multiplicative constant. This constant will be set by the behavior of the disorder 
on larger length scales. So if we want to find the value of a(x) at some point, we proceed as follows: 
consider short distance fluctuations to obtain a(x) as a solution of a wave equation on some short scale. 
The overall multiplicative constant on this wave function is not known. So, a(x) at a given point is 
known up to a multiplicative constant set by longer scales. Solving at a slightly longer length scale, we 
can find a(x) up to a constant set at even longer length scales. 

So, we would find that a(x) is given by a product of these multiplicative constants on larger 
and larger length scales, where these constants will be drawn from some random distribution. This is 
simply a statement that the zero energy wavefunction, or the value of the condensate, is a multiplicatively 
renormalizable operator, which gets a random multiplicative renormalization at each step. Thus, a(x) 
would have a log-normal distribution. Assuming scale invariance we would find that 
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a{x) = e p{x) (26) 
where, in order to produce the same fluctuations on all length scales, (3 has a Gaussian distribution with 

</?(?)/?(?)> = 4*(p + ?) (27) 
p a 

where g is some constant measuring the strength of the disorder. In two space dimensions, we can give 
(3 the probability distribution 



e 



i/d 2 *(9*/3(*)) 2 ( 28 ) 



Let us now derive this result from the RG. The zero energy wavefunction a(x) can be obtained 
by considering the correlation function of <j)(x, t) and <f>(y, t') where y is some point very far from x. 
So, we must look at how the correlation function G{x,t;y,t') renormalizes under the RG, beyond the 
simple calculation of the average correlation function in the previous sections. The interesting part is the 
renormalization of <j>(x). Let us look at this problem using moments. That is, we will look at the average 
over disorder of the n-th moment of G(x,t;y,t'). Let us start with the second moment for simplicity. 
This is an average of a product of two Green's functions. 

Consider what we may call a renormalization of a vertex. We must have both Green's functions 
starting at the same point, x,t. One may connect the two lines with a single scattering off of A, with 
momentum of order A running around the loop, with low momentum leaving the diagram. See figure 4 
for the appropriate diagram. The result of this is that the vertex is renormalized under RG flow as 

™-*L-^»°8 (29, 

where Vi is the vertex. 

In general, for the n-th moment, there are " ( -" 2 ~ 1 ^ ways to connect the lines at the vertex, and 
so we find 

d\vV n n(n-l) _ n(rc-l) 2 8 . 2C . . 

dmX " ~^— C3L ~ — 2~~ C3ClA S m 

Interpreting this result, we find that the n-th moment of a{x) at some point can be expressed, up to the 

multiplicative renormalization calculated here, as the n-th moment of a(x) at that point with a smoothed 

disorder potential. Looking at the renormalization, though, we see that it is exactly the form expected 

for a log-normal distribution of a(x). That is, if we took a(x) at some point to be expressed as a random 

multiplicative renormalization of the value of a(x) at that point calculated from a smoothed disorder 

potential, we would obtain a result of the form of equation (^0|) . 

Let us show that the moments are those expected from a log-normal distribution, as well as 

obtaining the value of g, measuring fluctuations in (3, as defined above. Using a(x) — e@( x \ computing 

the average of a n (x), over fluctuations in (3 down to some scale A, we find 

{a n (x))= [[d(3(p)}e^ dipinm ~£ m2) (31) 
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This is equal to 



e L dd ^ =e ^L d ^ (32) 

Finally, we see that 

-an (33) 



dlnA r(d/2)' 



So we find, comparing equations (30) and (p3|), that 



where we need to choose g to make the coefficients of the n 2 term in equations ( p0[ ) and (33) the same. 
The term in n can be different, as this simply represent an overall scale for the wavefunction. 

Given this form of the wave function, with log-normal fluctuations, the change in the low energy 
density of states should not be a surprise. Originally, the field <\> had the action 

d d x d {D ~ d) t (-d^x, t)a^0(ar, t) + A(a;)#) (35) 

Given the zero energy wave function, we can write this as 

d d xd D - d ha(x) 2 d„(a- x (x)^(x,t))d'"(a- 1 (x)(t>(x,t)) (36) 



This leads to a random stiffness problem, which is related to the problem of disordered SUSY quantum 
mechanics || and the problem of Dirac fermions in a random vector potential jj| . Both of these problems 
are known to have an increase in the low energy density of states. In particular, for the Dirac problem 
in which there exists a dimensionless measure of disorder, the low energy density of states is a power law 
with continuously variably exponent. 



VI. CONCLUSION 

In conclusion, we have presented a renormalization group treatment of the disordered large N 
sigma model. We have calculated the exponents describing the average correlation functions. However, 
there is still much that we would like to know. 

First, we have not considered the behavior of averages of higher moments of the correlation 
functions, or typical behavior of the correlation functions. These will be discussed in another publication 
H . The calculation of these higher moments is closely related to the calculation of the zero energy wave 
function. It is shown || that 

(G"(p, w)> oc (G(p, u))»p-»("-*»* L (37) 
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similar to the result for the wave function derived above. Here, C3 and L are the constants defined in 
section II, while the expectation value is an average of the n-th moment of the correlation function. 

Also, there exists another RG treatment of the same system g], based on a Callan-Symanzik type 
RG for the large N system. In that work, no perturbatively accessible fixed point was found, and I do 
not fully understand why their technique fails. Of course, this work [^| is not necessarily in contradiction 
with the results obtained here. Since no fixed point was found with their technique, this may simply be 
failure of technique, rather than indicating different results for the same system. 

Let us consider the differences in techniques used in more detail. It is conceivable that the lack of a 
perturbatively accessible fixed point in that work || is due to their use of a Callan-Symanzik technique, 
which is equivalent to the Wilson-Fisher RG only for a field theory which possesses a renormalizable 
continuum description. However, the 5-function interaction is non-renormalizable as a continuum theory 
(in that work a quartic interaction, which is renormalizable, was used, but the technique of summing 
polarization bubbles used there made the results equivalent to the ^-function theory) and the inverse 
polarization bubble used to calculate fluctuations in A from fluctuations in a 2 grows rapidly in the 
ultraviolet which can lead to poor behavior of the diagrammatic expansion of the continuum theory. A 
perturbative RG implies an ordering by momenta, in which high momentum processes are calculated 
before low momentum processes. For the expansion to be valid, this ordering by momenta must be 
correct, in that corrections to high momentum processes due to low momentum effects must be small. I 
have checked this for my procedure; I am not aware if it is true for other techniques on the same problem. 

Although results have been obtained for average quantities to all orders, one still may inquire 
about the radius of convergence of the expansion. It is apparent that if (D — d) > 2 then the dynamic 
critical exponent z resulting from equation ( |l8| ) becomes infinite. When (D — d) > 2, however, it is 
possible for the system to do a phase transition in a finite region in the d dimensional space, since there 
is enough volume in the remaining D — d dimensions. Thus, certainly when D — d = 2, and possibly 
before, the perturbation theory must break down. Further, results for the higher moments of correlation 
functions may have interesting behavior to higher order. 

Also, we would like to understand 1/N corrections to this problem. To lowest order in I/TV and 



e such corrections only modify the RG equations by changing the scaling dimension of a 2 1 1 1 in the pure 
system, and thus changing the equation to 

^£ = D - 2 + V - ( C3 - c 2 )c\K*- 2D S{l (38) 

= D - 2 + r, - (c 3 - c 2 ) Cl A 8 - 2 ^(l - B^*) - d/2 (39) 

where in the physically interesting case of d = 2, D = 3 we find r\ = -^z^- Then, we would expect to 
still find vd = 2 for the disordered critical point. 
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FIG. 1. Polarization bubble. Thick lines represent either scattering vertex off A or scattering vertex used to 
define a 2 in self-consistency equation. 




FIG. 2. Self-energy correction due to fluctuations in A. Joining the thick lines in a loop denotes averaging A 
over disorder in a 2 . Momentum of order A flows around loop. 




FIG. 3. Vertex correction due to fluctuations in A. This represents both renormalization of vertex defining 
scattering off of A and renormalization of vertex defining a 2 in self-consistency equation. 
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FIG. 4. Renormalization of vertex in computing higher moments of Green's function. Two Green's functions 
start at the same point. After Fourier transforming, this implies that they start with given total momentum. By 
including fluctuations in A, with momentum of order A running around the loop, one can define a renormalized 
vertex. 
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